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We obtain photometric constraints on the luminosity of the stellar remnant in SNR1987A using XMM-Newton and 
INTEGRAL data. The upper limit in the 2-10 keV band based on the XMM-Newton data is L x & 5 x 10 34 erg/s. 
We note, however, that the optical depth of the envelope is still high in the XMM-Newton band, therefore, this 
upper limit does not constrain the true unabsorbed luminosity of the central source. The optical depth is expected 
to be small in the hard X-ray band of the IBIS telescope aboard the INTEGRAL observatory, therefore it provides 
an unobscured look at the stellar remnant. We did not detect statistically significant emission from SN1987A in 
the 20-60 keV band with the upper limit of Lx 1.1 x 10 36 erg/s. We also obtained an upper limit on the mass 
of radioactive 44 Ti M( 44 Ti)< 10~ 3 M r 
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1. Introduction. 

Owing to its proximity, Supernova SN1987A in the Large 
Magellanic Cloud provides a unique possibility to study 
the physics of a supernovae and supernova remnants. 

Today it reveals itself as an extended object experienc- 
ing rapid evolution. The observed X-ray emission of the 
supernova is dominated by the interaction of the shock 
wave with the matter produced by the wind from a super- 
giant phase of its progenitor. As a result of this interaction 
the supernova was gradually increasing its luminosity dur- 
ing the past decade. This was followed by a strong bright- 
ening of the supernova on December 2002 whe n the shock 
wave e ncountered the so-called inner ring fsee lPark et alJ 
( 2004| and references therein) . 

Based on the Balmer lines in the optical spec- 
trum, i dentification of a mass ive blue supergiant pro- 
enitor jSonneborn et all Il987fl and the neutrino burst 
Hirata et all Il98*?fi coincident with the SN event, 
SN1987A was classified as a Type II core-collapse super- 
nova. The m ass of the progenitor was estimated to be 



15-20M r . 



[IT987lblSunvaev et a,lj|l98^ 

detected strong very hard X-ray source after about 169 
days from the explosion but this emission was connected 
with comptonisation of gamma-ray lines of Co 56 decay in 
the optically thic k envelope. The l i nes th emselves became 
visible later fe.g. ITeeearden et all Il989t) . 

The collapse of such a massive star should lead to the 
formation of a neutron star or a black hole. If a neutron 
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star was born with strong magnetic field and was rapidly 
rotating, it should appear as a luminous Crab-like X-ray 
source or as an accreting pulsar, which usually also have 
very hard X-ray spectrum with maximal contribution to 
luminosity in the INTEGRAL band. Neutron star with- 
out strong magnetic field or a black hole might become 
a bright X-ray source due to the a c cretio n of the matter 
of the dense envelope. IPark et all <)2004|) obtained from 
an imaging analysis of the Chandra data an upper limit 
Lx (2-10 keV)< 1.5 X 10 34 erg/s on the luminosity of the 
central source. However the envelope continues to be opti- 
cally thick and there is a possibility that its optical depth 
is still high in the standart X-ray band and does not per- 
mit us to observe the central source directly. It is impor- 
tant to check that it is not specially bright in the hard 
X-ray band, where the envelope is transparent both for 
the Thomson scattering and photoabsorption. In addition 
we decided to check the available data obtained by the 
XMM-Newton observatory, which has higher sensitivity 
for photons with energies ii 5keV than Chandra. 

Here we present results of analysis of INTEGRAL and 
XMM-Newton observations of SNR1987A and constrain 
the luminosity of the stellar remnant embedded in the 
envelope based on its X-ray spectrum. 

2. Data reduction and results 

2.1. XMM-Newton 

SNR1987A was observed by XMM-Newton in September 
and November 2000, April 2001 and May 2003. These 
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observations were processed with the Science Analysis 
System (SAS) v6.0.0. The data were screened for the soft 
proton flares by removing the time intervals where the 
count rate above 10 keV significantly exceeded the mean 
level. For the pointing in April 2001 because of numerous 
soft proton flares the background threshold was chosen 
higher than the quiescent background typical for time pe- 
riods without flaring activity. To be sure that this does 
not influences our results, we repeated analysis with dif- 
ferent thresholds on background level and a background 
subtraction from different regions on detector and found 
no significant variations in the obtained results. 

The angular resolution of the XMM-Newton, unlike 
that of the Chandra, is insufficient to resolve the remnant, 
and it appears as a point-like source on the extracted im- 
ages. However, XMM-Newton is more sensitive to the pho- 
tons with energy > 5 keV and therefore allows to recon- 
struct the sources spectra to higher energies. The EPIC 
spectra were extracted from a 30" — 40" circle regions 
around the SNR. Corresponding background spectra were 
extracted from the nearby regions. 

The spectra of the SN1987A remnant obtained with 
XMM-Newton on September 2000 and May 2003 to- 
gether with NEI model fit (the plane-parallel shock 
model l)Borkowski et alll2001af) where the plasma has not 
reached the collisional ionization equilibrium) are shown 
in Fig^ As is evident from the figure, the luminosity 
of the remnant increased by a factor of «3 during three 
years. This is expected from the interaction of the super- 
nova shock with the matter left from the stellar wind of a 
supergiant progenitor. It seems that there is some excess 
of hard X-ray emission in comparison with NEI model for 
pointing on May 2003. Such a high energy tail could be as- 
sociated with either a central source or a synchrotron radi- 
ation f rom the shell of the remnant (e.g. lBorkowski et all 
l2001b() . We focus on the possible association of the high 
energy excess with the central source. Detailed analy- 
sis of the supernova remnant spectra is not a subject o f 
this paper and cou ld be found in e.g. IPark et alJ (20041) : 
lAschenbach1l|2002l) . 

To obtain a conservative upper limit on the central 
source luminosity we shall assume that the total flux ob- 
served at the high energy part of the remnant's spectrum 
is produced by the central source. For this analysis we used 
the data of the early pointings (2000 and 2001 years), as 
a contribution of the emission originating from the inter- 
action of the shock with the medium would be smaller for 
them. Assuming the spectrum of the central source to be a 
powerlaw with a photon index 2.0, we obtained the upper 
limit (la) for its 2-10 keV luminosity L< 5 x 10 34 erg/s, 
assuming a distance to the source d = 50 kpc. The ob- 
tained value is in the same range, but somewhat higher 
than the upp er limit obtained by Chandra from the imag- 
ing analysis l|Park et alll2004h . 



2.2. INTEGRAL 

The internatio n al ga mma-ray laboratory INTEGRAL 
llWinkler et all Hp03) observed the Large Magellanic 
Cloud several times on January 2003 with a total ex- 
posure about 1 Mln sec. In our analysis we used pub- 
licly available data obtained by the ISGRI detector of 
the IBIS telescope of the observatory. ISGRI is sensi- 
tive to the photons in the E> 20 keV energy band. The 
data of INTEGRAL/IBIS /ISGRI were processed with the 
method developed by Eugene Churazov and described in 
Revnivtsev et al. (2004). Detailed analysis of Crab nebula 
observations suggests that with the approach and software 
employed, the conservative estimation of the uncertainty 
in measurements of absolute fluxes from the sources is 
about 10%. 

To investigate the hard X-ray emission from the 
SN1987A remnant we analyzed data obtained with 
INTEGRAL/IBIS /ISGRI in the 20-60 keV energy band, 
where the sensitivity of the ISGRI detector is maximal. We 
reconstructed the image of the central part of the Large 
Magellanic Cloud in this energy band (see Fig. EJ) and did 
not find statistically significant emission from the remnant 
with the upper limit of Fx 5s 3.7 x 10~ 12 erg/s/cm 2 (2a, 
assuming a Crab-like spectrum), that corresponds to the 
luminosity of Lx 5a 1.1 X 10 36 erg/s. 

To compare this result with results of Chandra and 
XMM-Newton observatories we can extrapolate it to the 
2-10 keV energy band. The spectra of pulsars in the su- 
pernova remnants are usually approximated by a power- 
law model with a slope of ~ 1.5 — 2. For such values of 
the photon index we can estimate the central source lumi- 
nosity in the 2-10 keV energy band corresponding to the 
INTEGRAL upper limit as L x 5S (0.6 - 1.6) x 10 36 erg/s. 

2. 3. Optical depth of the envelope 

The envelope of the supernova at early times was opti- 
cally thick both for photoabsorption and Thomson scat- 
tering. This, for example, follows from the observations of 
Rent gen observatory on th e KVA NT module of MIR sta- 
tion llSunvaev et aTlll987t see also lGrebenev fc Svunvaevl 
l|l987|) 1. As the envelope expands and its density de- 
creases, it s optical depth decreases wi th time as r oc t~ 2 . 
Given the lFransson fc Chevalier! ( 19871) model for the pho- 
toabsorption optical depth of the envelope, and assuming 
t~ 2 behaviour, we obtain T p hoto ~ 7 for a 5 keV photons 
at the epoch of XMM-Newton obervations. The photoab- 
sorption optical depth reaches unity only on energies ~11- 
15 keV. Thus, the optical depth is still high in the energy 
bands of Chandra and XMM-Newton. Therefore the up- 
per limits provided by these observatories correspond to 
the small fraction of the central source emission which 
escaped the envelope, while the real central source lumi- 
nosity could be much higher. 

Similar estimate of the optical depth for 40 keV pho- 
tons give value Tphoto ~ 0.02, therefore the hard X-ray ra- 
diation escapes the envelope without absorption. For the 
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Fig. 1. Left: EPIC PN spectra of the SNR1987A. The upper curve corresponds to pointing on May 2003, the lower 
to pointing on September 2000. Right: Deconvolved photon spectrum of SNR1987A on April 2001 with the powerlaw 
model describing the data on energies above 4keV. 



Thomson optical depth we obtained also very low values, 
tt ~ 0.005 — 0.1 depending on the model of matter distri- 
bution in the envelope. Was the tt close to unity or even a 
bit higher, it would not affect the flux in 20-60 keV energy 
band because of small e nergy change for one sca ttering, 
AE/E » hv/m e c 2 ~ 0.1 i|Pozdniakov et al.lll983|) . Thus, 
the INTEGRAL observations offer us a direct look at the 
central source and its upper limit constraints the real lu- 
minosity of the stellar remnant. 

Absorbed emission in the standart X-ray band is 
reemitted in the optical and infrared bands and an upper 
limit on the central source luminosity could be estimated 
based on a bolometric lum inosity of the SN R in these en- 
ergy bands. For example, ISuntzefj l|l997l) estimated the 
bolometric luminosity of the remnant log(L)~ 36.1 — 36.4 
after the 3600 days from the explosion. This value could be 
interpreted as an upper limit on the central source lumi- 
nosity corrected for the absorption, what is close to a limit 
obtained from the INTEGRAL data. Obviously, the up- 
coming Spitzer data will provide an interesting estimates 
on the luminosity of the central source, absorbed and rera- 
diated in the infrared band by the optically thick for the 
standart X-rays envelope. 

2. 4- Limit on the nucleosynthesis of M Ti 

While the early (except for the first few weeks) bolo- 
metric luminosity of the supernova was powered by the 
56 Co decay (ti/ 2 ~77 days), at later times elements with 
longer decay times, like 44 Ti (ti/ 2 ^60 years), give the 
main contribution. Decay of the 44 Ti gives rise to two 



gamma photons with energies 67.9 keV and 78.4 keV. 
Simple estimations on the fl uxes of these lines give value ~ 
3 • 10- 13 erg/s/cm 2 (see e.g. iMotizuki fc Kumaeal l2004|b 
assuming mass of titanium M( 44 Ti)= 10~ 4 Mq. To esti- 
mate the mass of the synthesed 44 Ti, we calculated fluxes 
observed by INTEGRAL in the 61-73 keV and 73-88 keV 
energy bands, F x (61-73 keV)=4.30±1.25-10- 12 erg/s/cm 2 
and F x (73-88 keV)=2.50 ± 1.60 • 10- 12 erg/s/cm 2 . In the 
first case the flux has a formal significance of 3.4<r, how- 
ever we consider it insufficient to claim a reliable detec- 
tion of the line emission. Therefore to estimate the mass of 
44 Ti, we assume both fluxes to be upper limits and obtain 
M( 44 Ti)^ 1O~ 3 M . 

3. Conclusion 

Based on the data of XMM-Newton and INTEGRAL, we 
provide broadband spectroscopic constraints on the lumi- 
nosity of the central source in the SN1987A. Assuming 
that a hard energy tail of the XMM-Newton spectrum 
of the remnant is dominated by a central source, we ob- 
tained an upper limit L< 5 x 10 34 erg/s for 2-10 keV cen- 
tral source luminosity uncorrected for absorption, which 
is somewhat higher than the upper limit obtained by 
Chandra from imaging analysis. INTEGRAL did not de- 
tect any statistically significant emission from the remnant 
and provided an upper limit for 20-60 keV luminosity of 
Lx ^ 1.1 x 10 36 erg/s, what corresponds to 2-10 keV lumi- 
nosity Lx ^ (0.6 — 1.6) x 10 36 erg/s assuming Crab-pulsar 
like spectrum. The obtained upper limits on the possible 
central source luminosity in comparison with the typical 
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Fig. 2. Map of the part of LMC around of the SN1987A position, obtained with INTEGRAL in the 20-60 keV energy 
band. Position of SN1987A is shown by the cross. 



spectra of Crab-like pulsar, accreting X-ray pulsar and 
accreting black hole are shown on Fig. [3] 

Based on the existing calculations of the photoabsorp- 
tion optical depth of the envelope, we show that absorp- 
tion in the standart X-ray band of the XMM-Newton and 
Chandra is still high and obtained upper limit corresponds 
to luminosity not corrected for absorption, which could 
be only a small fraction of a real luminosity. The opti- 
cal depth in the INTEGRAL energy band is significantly 
less than unity and obtained upper limit constrain real 
luminosity of a stellar remnant. 

We also obtained an upper limit on the mass of 44 Ti 
M( 44 Ti)^ 10^ 3 M Q . 
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Fig. 3. The typical spectra produced by accreting X-ray pulsar, accreting black hole l|Gilfanov et all l2000h and 
rotation-powered Crab- like pulsar. All spectra are normalised to the lumin osity of 2.5 • 10 36 erg/s in the 1-12 keV energy 
band. This corresponds to the bolometric optical and infrared luminosity l|Suntzefj . Il997|) , which could be interpreted 
as an upper limit on the central source luminosity absorbed by the envelope. Upper limits on the luminosity of central 
source in SNR1987A obtained by the XMM-Newton and INTEGRAL observatories are also shown. 
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